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ABSTRACT 

Pointed  bubbles  with  a  120°  angle  at  the  apex  rising  in  a  two-dimensional 
tube  are  considered.  A  solution  is  found  at  a  unique  Froude  number  of  0.36. 
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POINTED  BUBBLES  RISING  IN  A  TWO-DIMENSIONAL  TUBE 


Jean-Marc  Vanden-Broeck 

We  consider  a  bubble  rising  at  a  constant  velocity  U  in  a  two- 

dimensional  tube  of  width  h.  We  choose  a  frame  of  reference  moving  with  the 

bubble  and  we  assume  that  the  bubble  extends  downwards  without  limit.  We 

denote  by  2y  the  angle  between  the  free  streamlines  at  the  apex  of  the 

bubble  (see  Figure  1).  We  shall  neglect  the  effect  of  surface  tension. 

Therefore  the  problem  is  characterized  by  the  Froude  number 

F  =■  U/(gh)1/2  .  (1) 

Here  g  denotes  the  acceleration  of  gravity.  This  problem  was  considered 

before  by  Garabedian 1 ' * ,  Vanden-BroeckJ  and  Modi*. 

*  v 

Garabedian1  presented  analytical  evidence  that  a  solution  with  y  ■  ^ 
exists  for  each  value  of  the  Froude  number  F  smaller  than  a  critical  value 
F  •  Vanden-Broeck  solved  the  problem  numerically  and  found  F_  =  0.36.  In 

C  C 

o 

addition  Vanden-Broeck  showed  numerically  that  a  solution  with  y  “  0  (i.e. 

with  a  cusp  at  the  apex)  exists  for  each  value  of  F  greater  than  Fc> 

2  4  IT 

Garabedian  and  Modi  pointed  out  that  solutions  with  y  *  —  might  also 

exist.  The  flow  in  the  neighborhood  of  the  apex  would  then  be  similar  to  the 

classical  Stokes  flow  near  the  crest  of  the  highest  progressive  wave.  Other 

flows  with  an  angle  of  120°  on  the  free  surface  have  been  found  by  Vanden- 

Broeck  and  Keller5  and  by  Goh  and  Tuck6. 

In  this  paper  we  present  numerical  evidence  that  a  solution  with  y  -  y 

exists  only  for  F  *  Fc  ~  0.36. 
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We  define  dimensionless  variables  by  choosing  h  as  the  unit  length 
and  U  as  the  unit  velocity.  We  introduce  the  potential  function  e  end  the 
stream  function  Without  loss  of  generality  we  choose  <p  =  0  at 

x  ■  y  ■  0  and  *  0  on  the  bubble  surface.  We  denote  the  complex  velocity 

by  c  -  u  -  iv  and  we  define  the  function  T  -  i0  by 

C  -  u  -  iv  «  eT”i8  .  (2) 

7  1 

Following  Birkhoff  and  Carter  and  Vanden-Broeck  ,  we  define  the  new 
variable  t  by  the  relation 

e"*f  -  ~  (t  +  .  (3) 

Here  f  =  <p  +  i<J>  denotes  the  complex  potential.  It  can  be  shown  that  this 
transformation  maps  the  flow  domain  onto  the  unit  circle  in  the  complex  t- 
plane  so  that  the  free  surface  SJ  goes  into  the  portion  of  the  circumference 
in  the  first  quadrant  (see  reference  3  for  details).  The  complex  t-plane  is 
shown  in  Figure  2. 

The  flow  in  the  neighborhood  of  the  apex  S  of  the  bubble  is  locally  a 
flow  inside  an  angle  it  -  y  (see  Figure  1).  Therefore  Z  vanishes  like 
fY/,1T  as  f  0.  Using  (3)  we  have 

Z  ~  (1  -  t)2Y/ir  as  t  -*•  1  .  (4) 

Following  the  procedure  outlined  in  reference  3  we  represent  Z  by  the 
expansion 

Z  *  -[-in  C( 1  +  t2)]1/3(-in  C)"1/3(1  -  t2)2r/w[1  +  l  ant2n]  .  (5) 

n=1 

Here  C  is  an  arbitrary  constant  between  0  and  0.5.  We  chose  C  =  0.2.  The 
coefficients  ^  have  to  be  determined  to  satisfy  the  condition  that  the 
pressure  is  constant  on  the  free  surface  SJ.  We  use  the  notation  t  =  |t|e*a 
so  that  points  on  SJ  are  given  by  t  =  e*CT,  0  <  a  <  -j*  Vanden-Broeck3 


showed  that  the  pressure  condition  on  SJ  can  be  written  as 


„  2t  dr  1  -t  2 
it  cot  oe  —  +  — r  e  cos  6  =  0. 
do  _2 
F 


(6) 


Here  t(o)  and  6(a)  denote  the  values  of  t  and  6  on  the  free  surface 
SJ. 

For  y  =  0  and  Y  =  "J/  (5)  is  equivalent  to  the  representations  (15), 

(18)  and  (19)  of  reference  3. 

w 

In  order  to  find  solutions  with  a  120#  angle  at  the  apex,  we  set  y  = 
in  (5)  and  determine  the  coefficients  an  by  using  the  procedure  outlined  in 
reference  3.  Thus  we  truncate  the  infinite  series  in  (5)  after  N  -  1 
terms.  We  find  the  N  -  1  coefficients  an  and  the  Froude  number  F  by 
collocation.  Therefore  we  introduce  the  N  mesh  points  0j  *  (I  -  “) , 

I  =  1 , . • « ,N  and  satisfy  (5)  at  the  N  points  a^.  T^iis  leads  to  a  system 
of  N  nonlinear  algebraic  equations  for  the  N  unknowns  aR,  n  *  1,...,N  -  1 
and  F.  We  solve  this  system  by  Newton's  method.  Once  this  system  is  solved, 
the  shape  of  the  bubble  is  obtained  by  numerically  integrating  the  identity 

(7) 

In  Table  I  we  present  numerical  values  of  F  corresponding  to  different 
values  of  N.  The  value  0.3577  appear  to  be  correct  to  4  decimal  places.  The 
corresponding  profile  of  the  bubble  is  shown  in  Figure  3. 

As  a  check  on  our  results,  we  modified  the  scheme  to  allow  the  Froude 
number  F  to  be  prescribed.  Solutions  could  not  be  obtained  for  other  values 


1*  +  4  lX  =  _-T+i8 
3<p  1  3<p  e 


N 


10  0.35733 

30  0.35763 

50  0.35771 

70  0.35775 


Table  1:  Values  of  the  Froude  number  F  corresponding  to  the  bubble  with  a 
120°  degree  angle  at  the  apex  for  various  values  of  N. 


The  results  in  reference  3  results  and  the  present  calculations  show  that 


TV 


y  -  t 


2  ' 
TT 


0  <  F  <  F  , 
C 


y  =  3  '  F  -  Fc  , 


Y  =  0, 


Fc  <  F  <  ob 


(8) 

(9) 

(10) 


where  F  ~  0.36. 

C 

2 

The  results  (8)— (10)  do  not  support  Garabedian's  conjecture  that  bubbles 
with  a  120°  angle  at  the  apex  exist  for  all  F  between  0.23  and  0.36. 

D 

Finally  let  us  mention  that  Vanden-Broeck  solved  the  problem  with 
surface  tension.  He  found  that  for  each  value  of  the  surface  tension  there 
exists  a  countably  infinite  number  of  solutions  for  which  y  =  -j  .  As  the 

surface  tension  tends  to  zero,  all  these  solutions  approach  a  unique  limiting 

* 

configuration  characterized  by  F  =  F  ~  0.23.  The  profile  corresponding  to 

* 

F  =  F  was  found  to  be  in  good  agreement  with  experimental  data. 
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